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Abstract 

A model consisting of quintessence scalar field interacting with cold 
dark matter is considered. Conditions required to reach Wd — are 
discussed. It is shown that depending on the potential considered for 
the quintessence, reaching the phantom divide line puts some con- 
straints on the interaction between dark energy and dark matter. This 
also may determine the ratio of dark matter to dark energy density at 
Wd = -1. 

PACS: 98.80.-k, 95.36. -l-x 



1 Introduction 

Dark energy model is a candidate to explain the present accelerated expan- 
sion of the universe [I]. In this scenario 70% of the universe is assumed 
to be permeated by a smooth energy component with negative pressure, 
dubbed as dark energy. A simple model introduced to describe dark energy 
is a scalar dynamical field with a suitable potential [2]. In this model an 
exact solution for Friedmann equations is only accessible for some kinds of 
potentials [3]. 

Based on recent astrophysical data, it seems that the dark energy com- 
ponent has an equation of state parameter, Wd — —1, in the present epoch 
[1]. Therefore to study the dynamical behavior of dark energy in the present 
era, instead of trying to find an exact solution, one can restrict himself to 
the region Wd — —1, where equations can be solved approximately. 

In a noninteracting dark energy model, it is expected that density of 
dark matter (pm), in the present epoch, be very less than the density of 
dark energy component (pd)- This lies on the fact that the equation of state 
(EoS) parameter of dark energy, Wd, is less than — therefore it redshifts 
more slowly than the (dark) matter component. But the ratio of dark matter 
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to dark energy, r, is of order unity in the present epoch: r ~ 3/7, this is 
known as the coincidence problem [5] . This forces us to consider interaction 
between dark energy and matter which ahows energy exchange between 
these components [6j. 

In this paper we assume that the universe is filled with a scalar field, 
whose EoS satisfies —l<Wd<—^ (dubbed as quintessence), and (cold) 
dark matter with mutual non- gravitational interaction. As we have men- 
tioned obtaining an exact solution for the Friedman equations in the pres- 
ence of interactions even for simple potentials is not straightforward. Hence 
using some natural conditions which must be obeyed by the quintessence 
field in the vicinity oi Wd — —1, we seek the required conditions (posed 
on the interactions and scalar potentials) allowing the system to reach at 
Wd = —1 without requirement to obtain an analytic and an exact solution 
for the problem. So using this method one can show that if a proposed 
quintessence model with a specific potential and interaction permits reach- 
ing to ~ — 1 in the present era or not. 

Besides, conditions on the interaction term and the potential which may 
be functions of energy densities may pose some natural conditions on the 
ratio of dark matter to dark energy density showing whether the occurrence 
of coincidence problem near the era where t/;^ = — 1 is occurred is only a 
coincidence. 

At the end, using the interaction Q = H{XmPm + ^dPd) and via Taylor 
series expansion we illustrate and emphasize our results. 
We use units G = ks = c = 1 throughout the paper. 



2 Quintessence model and Wd = 

The spatially flat Friedmann Robertson walker space time is described by 
the metric 

ds^ = df - a^{t){dx^ + dy"^ + dz^), (1) 

where a{t) is the scale factor. We assume that this universe is filled by dark 
energy and (cold) dark matter with densities pd and pm respectively. Dark 
energy component is a scalar field, (/>, with potential 1^ (</))• Energy density, 
> 0, and pressure, < 0, of dark energy are given by 

Pd = \ + y{<t>) 

Pd = f (2) 

The EoS parameter of dark energy, given by 

Pd 4>^ - 2V{<j,) 
Wd = — = -■ , (3) 
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satisfies —1 < Wd < —5. The scalar field with this EoS parameter is dubbed 
quintessence. The inequality < is necessary for accelerating the ex- 
pansion of the universe. The Hubble parameter, H = ^, satisfies Friedmann 
equations 

3 

H = -MP + p), (4) 

where p and P are the total energy density and pressure of the universe: 
p = P(i + Prm P = Pd- We consider an interaction between dark matter and 
dark energy components: 

Pd + 3H{Pd + Pd) = -Q 

p'm + 3Hpm = Q. (5) 

By substituting in the above equation we find the evolution equation for 
the quintessence field 

<P{4> + 3H<j) + V'{^)) = -Q. (6) 
The ratio of dark matter to dark energy, r = — , satisfies 

r = 3Hrwd+ —{1 + r). (7) 
Pd 

Using Wd < and the above equation one can show that in noninteracting 
quintessence model: r < —Hr, and therefore r — )• eventually, in contrast 
to the recent data which assess r ~ 0{1). We can also obtain an expression 
for time evolution of r in terms of fid defined by = ^ (where pc = 
is the critical density): 

By comparing ([7]) and ([8]) we arrive at 

= ^ 9. (9) 

SHQdil-^d) SHpdil-ndY 

In the absence of interaction we have ild > H{^~^d)y which implies that Qd 
is an increasing function of time. In contrast, flm, defined by = is a 
decreasing function in this case. Hence the introduction of Q in the above 
equation, via exchanging energy, depending on the form of Q may prevent 



r 



%^ to go to zero, and can be a remedy for the coincidence problem. 



Eq. ([9]) is a general equation for dark energy models. In the scalar field 
quintessence model we can also make use of 

2V{cP) = il-Wd)pd 

= {1+Wd)pd, (10) 
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to obtain 

[(1 - Wd)pd] = 2V'{<I^)^, (11) 

and subsequently 

± ^(1 + Wd)pdV'{4') = \{-WdPd + (1 - Wd)pd). (12) 
+ (— ) corresponds to > (<)0. By substituting 



Pd 



Svr 



^H^dil + WdQd) + ^d) , (13) 



which can be verified by taking the time derivative of pd = ^dPc into ()12p . 
and by making use of ([9]), we arrive at 



Wd = t2V'{4>).^^ - 3H{1 - wj) - -^(1 - Wd). (14) 

V Pd Pd 

We use the above equation to study the behavior of the model in the vicinity 
of the time, t = to, when Wd{tQ) = —1. The equation of state parameter 
of the quintessence is equal to or greater than —1, Wd > —1, therefore at 
to we must have Wdito) = 0. Otherwise there will be a neighborhood of t^, 
in which Wd < —1- Hence a necessary condition to reach at Wdito) = — 1 
in the quintessence model is = at t = to- Even when to ^ oo, this 
is asymptotically valid. This lies on the fact if — )• — 1 when t — )• cxd, we 
must also have ^(g 5R) = 0. For lim^^_^_i -^^^{1 + Wd) = (e.g. for 

bounded ^J^-^ at t = to) (HH) reduces to 

—— = 0. 15 

Pd{to) 

E.g. if one takes the interaction term such as Q = )^PmPd ( A > 0), as in 
the present epoch r ~ 0(1), we have ^ 7^ and it is clear that Wd = — 1 
cannot occur in the present era. 

Note that because of the presence of additional terms in ([6])with probable 
singular behavior at the limit Wd — >• —1, Q{to) = (at Wd = —1) may not 
be derived directly from ([6]) unless lim^^_j._i ^J^^ \/(l + '^d) = 0. 

Beside the above natural condition, Wdito) > is also a necessary con- 
dition to reach at Wd{to) = —1- This is due to the fact that Wd = —1 may 
only be the global minimum of Wd{t). Note that even if Wd = —1 occurs at 
to — >• oo, limt^ooWdi't) = must be satisfied (provided Wd G JR). Now we 
must explain this natural condition in terms of the interaction term and the 
potential (inconsistency of the potential and the considered interaction with 
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this condition indicates that the model is not able to reach at Wd = —1)- To 
do so by getting another time derivative of eq. (jl4p . at t = to we obtain 

= -2^^ T 2 . (16) 

Note that the right hand side must be evaluated at t = to- By putting dH]) 
into (fT6]) we get 



Wdito) = -2 +2 ± 2 17 



But at t = to 



Pd 



— + — 



dt Pd 
Q 
Pd 

holds, therefore Wd(to) > requires that at t = to: 



(18) 



-Q{to) + {V'm^±^-^>0. (19) 
V^F(0) Vl + Wd 

For interactions of the form Q{pm, Pd)^ at t = to we have <5(to) = —3Hpm-§^ 
and the equation at t = to reduces to 

-,Hp^^HV'm'± ^Jfi /i >0. (20) 

OPm \/^((A) V 1 + 

Note that the form used in ()17p for the third term in the right hand 
side is suitable only when Q can be expressed in terms of 1 + Wd, (e.g. 
for interactions containing positive power of (p which following (jlOp can be 
written in terms 1 + Wd like the interaction considered for the inflaton 
during the reheating process :Q cc Pd + Pd = 9^)- In general, at t = to, using 
I'Hopital's rule, it is also possible to write (jl7p in the form 



y'(0) 



Wd{to) = -2^^-j^^V^)- (21) 

For intermediate time (we mean to ^ oo), and for t < tg, we have tf)^ < and 
for t > to, > holds. {+(—)} corresponds to the case where (j) < (^)O 
when t < tf). ()2ip has real roots (for liid) provided that at t = tg 

4Q(to)<(y'(0))'. (22) 
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The above inequality can be viewed as a constraint on the parameters of the 
model. In general, if the first nonzero derivative of Wd at to is of order n, at 
t = to we have 




In this case, evenness of n and wj^ (to) > 0, together with Wd{to) = are suf- 
ficient conditions for Wd to have a global minimum at to- The generalization 
of (HI]) is then 

^''^ - UJ ^ '(f^ It^J V ^''^ 

To derive the above equation we have assumed that ^ ^J^ ^ and its time 

derivatives up to order (^ — 1) are continuous and bounded at t = to- If 
Wd tends asymptotically to —1, all of derivatives of Wd may be zero in this 
limit. In this situation to is the point of infinite flatness and for inflnitely 
differentiable Wd, can only occur at infinity, to — s- oo. 

In the following, to elucidate our results, as an example, we consider the 
interaction [7J 

Q = H{\mPm + ^dPd)- (25) 
Following (|15p . we deduce that in order to reach at = — 1, we must have 

r(to) = (26) 

This equation determines the ratio of dark matter to dark energy in terms 
of the parameter of interaction at t = to . As we consider r as continuous 
function of comoving time, we also expect that this ratio is approximately 
given by r ~ — ^ in the vicinity of to where Wd{tQ) = —1 occurs. 

As a result, near Wd = —1, the value of r(to) is specified by the constant 
parameters of the interaction. If the present value of Wd is believed to be 
it)^ PS — 1 [4J, based on astrophysical data [5], we can get — ^ ~ 3/7. In this 
method we cannot assess and \m separately. Note that for interactions 
which do not satisfy ([26]) (e.g. models with ^ > 0), tt^d = — 1 is not 
accessible. 

To investigate the condition ()22p . we note that at tf^ = — 1 

g(to) = 8^Ad(l + r)p2. (27) 

Hence the inequality (j22]) becomes 

d,2n\d{l + r)pi{t^)<{V'm\ (28) 
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Note that the above constraint depends on the form of the potential of the 
quintessence field. We shall discuss this issue shortly after some remarks. 

The conditions like (jlSp and ()22p . are only necessary conditions to reach 
Wd = —1. Indeed we have used the fact that if in an interacting quintessence 
model 11!^, = —1 is achieved, then equations like [15] and [22] must be hold. But 
we did not prove that Wd = —1 is allowed in the model. In fact describ- 
ing the exact form of Wd (to see whether Wd = —1 is accessible) requires 
solving the equation ([6|) with one of the Friedmann equations in which 
in the presence of interaction (even in its absence), as we mentioned in the 
introduction, is not straightforward for a general potential. Despite this, 
to study more about the behavior of the system near Wd = —1, we can re- 
strict ourselves to the neighborhood of t = to where the Hubble parameter 
is presumed to be differentiable and consider the interaction (j25p . For a 
differentiable (at least in an open interval containing t = to) Hubble param- 
eter, in the vicinity of to (an open set containing to) following [8] we use the 
following Taylor expansion 

H = ho + hiit - tof + Oit - tof^\ 13 > I. (29) 
In this open set we consider the following expressions (Taylor expansion): 

Wd = -i+wo{t-tor + o{t-tor+^ 

r = ro + nit-toT' + 0{t-tor+\ (30) 

The equation of state parameter of the system is related to Wd through the 
relation w = Wd^d- We have also w = —1 — j^^, therefore near Wd = — 1, 
w = —0,d and / 0, and /3 in ([29]) begins with (3 = 1. We write the 
equation ([7]) as 



3rH 



Wd + U^^] (Xd + rXm) (31) 



Putting ([29]) and ([30]) into ([3T]) gives: 7 = 1 and ri = — 3ro/io. To elucidate 
our results we must specify the potential. Here we consider the quadratic 
and the exponential potentials. For the quadratic potential 

V{(t>) = \m^4>\ (32) 
by substituting ([29]) and ([30]) in ([Tl]) we arrive at 



± 2m./(l - vol) = Wd + 3H{1 - wj) + H{Xd + A„r)(l - Wd). (33) 



Comparing the coefficients of the expressions with the same power of t in 
both sides of ([33]) forces us to take 
([T3]) . We also obtain the equation 



both sides of ([33]) forces us to take a = 2 and ro = — , in accordance with 



± 2m\/2wo = 2wo + 2\mriho 

= 2wo-Q\dhl (34) 
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in agrement with the previous result ()2ip . The above equation has real roots 
provided that 

8^Ad(l - ^)<P\to) < 1, (35) 



which is the same as (|28p . This poses a condition on the value of the 
quintessence field at Wd = —1- In terms of total energy density this in- 
equality may be written as 

167rAdp(to) < (36) 
For the exponential potential 

y = 7;oexp(A(/)), vq > 0, (37) 

(fHl) reduces to 



T A(l - Wd)Vpdil + Wd) = Wd + 3H{1 - wj) + H{Xd + Amr)(l - Wd). (38) 

Again, by substituting ([2^]) and ([5n|) in ([55]) . and by comparing the coeffi- 
cient 
and 



cients of the same power of t in both sides we arrive at: a = 2, vq = —j^ 



=F ><VPd{'to)wo = wo + 3Xdhl, (39) 



in agrement with (|15p and (|2ip . The necessary condition to have real roots 
for (p9]l is then 

327rArf(l + ro) < A^ (40) 
By taking rg = —j^ ~ 3/7 as the present estimated value, we obtain 

3207r . . 2 
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Ad < A^ (41) 



which shows that access to it;^ = — 1 is not feasible for models with arbitrary 
A and Xd- 

3 Conclusion 

In this paper we considered a spatially flat FRW universe composed of dark 
matter and dark energy components. The dark energy was assumed to be 
a quintessence scalar field interacting with dark matter (see ([2|) and ([5])). A 
general expression for time derivative of EoS parameter of dark energy was 
derived (see (dH)), upon which we discussed some necessary conditions and 
relation between the interaction term and the potential of the quintessence 
to reach at Wd = —1 (see (fT5]) and ([22]) ). We also examine our results 
by approximation method based on series expansion near the time when 
Wd = -1. 
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